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ABSTRACT 

Theoretical models predict that some of the first stars ended their lives as extremely energetic pair- 
instability supernovae (PISNe). With energies approaching 10'^^ ergs, these supernovae are expected 
to be within the detection limits of the upcoming James Webb Space Telescope (JWST), allowing 
observational constraints to be placed on the properties of the first stars. We estimate the source 
density of PISNe using a semi-analytic halo mass function based approach, accounting for the effects 
of feedback from star formation on the PISN rate using cosmological simulations. We estimate an 
upper limit of ~0.2 PISNe per JWST field of view at any given time. Feedback can reduce this 
rate significantly, e.g., lowering it to as little as one PISN per 4000 JWST fields of view for the 
most pessimistic explosion models. We also find that the main obstacle to observing PISNe from the 
first stars is their scarcity, not their faintness; exposures longer than a few times 10^ s will do little 
to increase the number of PISNe found. Given this we suggest a mosaic style search strategy for 
detecting PISNe from the first stars. Even rather high redshift PISNe are unlikely to be missed by 
moderate exposures, and a large number of pointings will be required to ensure a detection. 



1. INTRODUCTION 

Understanding the formation of the first stars and 
galaxies is one of the central challenges of modern cos- 
mology, as they mark a significant increase in complex- 
ity from the simple initial conditions of t he pri mordial 
universe during the 'dark ages' (e.g., [Barka na fc Loeb| 
200T1 IMiralda-iSsclIdgl [20031 iBromm et anT l2009; ToeU 
2010 1 . The basic properties of these so-called Popula- 
tion III (Pop III) stars have been reasonably well estab- 
lished, with the consensus that the first stars formed in 
dark matter ' minihalos' on the order of 10^-10^ Mf?) at 
high redshifts (Couchman & Rees 1986 
[19961 



Numerica 



Haiman et al 
simulations ot 



Tegmark ct al. 1997) 
fhc collapse of primordial metal- free gas into these halos, 
where molecular hydrogen is the only available coolant, 
had suggested that the first stars were predominantly 
very massive, with M* > 100 Mf?) and a to p-heavy i nitial 
mass function (IMF) (e.g.,|Bromm et al.|1999 , 20 02l|Abel| 
et aL||2002 '; Bromm fc Larson||20U4| | Yoshida et al.f|2006[ 
O'Shea fc Norman, i2007D. More recent work has found 



significant fragmentation (jStacy et al.||2010 Clark et al. 


2011 Greif et al. 2011 


2U12p, and experienced strong 


protostellar feedback (F 


iosokawa et al.||2011 Stacy et al. 



star formation, with lower characteristic masses (on the 
order of 50 rather than 100 M©) and a much broader 
IMF now expected. The light from these stars ended the 
dark ages and fundamentally transfo rmed the univ erse, 
beginning both the reionization (e.g, Meiksin||2009 ) and 
the chemic al enrichment of the universe (e.g, Karlsson 
eFaLllmTl ) 



Given the top-heavy nature of Pop HI star formation 
and the fact that low-metallicity stars are unlikely to un- 
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dergo significant radiatively driven mass loss (Kudritzki 



2002 1 , one interesting possibility is that some of the first 
stars died as pair-instability supernovae (PISNe), a sce- 
nario that has significant consequenc es for the chemi- 
cal enrichment history of the universe (|Heger fc Woosley 



2002|[Tumlinson et al. 2004; Karlsson et al.|2008| ). Basic 



one-dimensional models predict that stars with masses 
in the range 140-260 Mq will under go a pair-production 



instability a nd explode completely (Barkat et al. 1967 



Fraley 1968) . During core oxygen burning, a combination 
ot high temperatures and relatively low densities results 
in the formation of pairs, removing pressure support 
from the core. Following the subsequent contraction, the 
ignition of explosive oxygen burning completely disrupts 
the progenitor, resulting in a significant contribution to 
the metal enrichment of the surrounding medium. More 



recently, Chatzopoulos fc Wheeler (2012) and Yoon et al 
(20121 have found that stars with initial masses as low 



as 65 JVIq can encounter the pair-production instability if 
they are rapidly rotating. In this scenario, strong rota- 
tionally induced mixing causes nearly homogeneous evo- 
lution, such that the star is converted almost entirely 
to helium before the next phase in its evolution begins. 
These extremely energetic explosions — approaching 10^'^ 
ergs for the most massive models — are very luminous, in 
part due to the large amount of ^^Ni produced, and are 
also very tem porally exte nded as a result of the large 
mass ej ected (Fryer ct al. 2001; Heger & Woosley 20021 
Heger et al.,2003j Joggerst fc Whaien,|2011, ,Kasen et all 
201l| : ' 



Xsecond possiblility for reaching such extreme explo- 
sion energies in slightly lower mass stars is the hypernova 
scenar io for rapidly rotating star s that undergo core col - 
lapse (Umeda fc Nomoto 2003, |Tominaga et al. 2007). 
During the collapse, accretion onto a central black hole 
powers a jet which induces a highly energetic explosion. 
While recent work has decreased the expected mass of 
the first stars, they have also bee n found to rotate m ore 
rapidly than previously thought (Stacy et al. 2011), in- 
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creasing the plausibility of this scenario. Evidence sup- 
porting this hypothesi s has recently been presented by 
Chiappini et al.| ( |2011[ ). 

While not an example of a Pop III star, the recent dis- 
covery of the extremely luminous supernova (SN) 2007bi, 
identified as a possible PISN, i n a metal- poor dwarf 
galaxy at a rcdshift of z ~ 0.1 (Gal- Yam et al. 20091 
suggests that PISNe may be possible in the local uni- 
verse u nder rare circumstanc es. Further supporting this 
picture, [mosleyet al.| ( |2007[ ) have shown that SN 2006gy 
( Smith et al.|2U07 ) is well modeled by a pulsational pair- 



the range of theoretical PISN models using the i mplicit 



instability model. For stars with initial masses in the 
range ~1OO-14OM0, the star encounters the pro- 
duction instability, but the resulting explosive ignition 
of oxygen is insufhcient to unbind the star. Instead it 
ejects a shell of material before settling back into a sta- 
ble configuration. The star encounters this instability 
several times until the mass of the helium core drops be- 
low ^40 M0 , after which the star can proceed to silicon 
burning and eventually undergo core-collapse. 

With the upcoming launch of the James Webb Space 
Telescope (JWST) we will be able to probe the epoch of 
first light in unprecedented detail. While the first stars 



themselves are unlikely to be visible (e.g, Bromm et al 
2001| [Pawlik et al.|201ip , some of the SNe that end their 



lives should be with in the detection limits of the JWST 
'e.g., Mackey et al. 



2003 Scannapieco et al.|2005 Gard 



ner et al. 2006). While the basic properties of PlSJNe and 
the effec t they have on their environment has bee n well 
studied (|Mori et al. 12002; Bromm ct al. 2003, Furlanetto 



fc LoeblpOOl^l IR itayama & Yoshida 2005 jWfia^eFaT 
20081 IWise & Abeiii2008( iGreif et al.ii2010l), the source 



Abeln2008( |Greif et al.||2010^ 
density of these events has yet to be well constrained. 

The first attempt at estimating the number and ob- 
servability of SNe at high rcdshift was made by |Miralda-| 
Escude & Rees ( 1997 ) , who calculated the all-sky SN rate 
based on estimates of the total metals produced by a typ- 
ical SN and the observed metallicity of the intergalactic 
medium (IGM) at high redshifts. This yields ~1 SN 
yr"""^ arcmin"^ ai z ^ b. Other early work attempted to 
model the SN rate based on the empirically determined 
star formation rate out to high redshifts (z ~ 5; e.g., 
Madau et aL]|1998[ [Pahlen fc Fransson||1999[ ) . It should 
be noted that these attempts focused on i'ype II SNe, not 
PISNe, whic h are the focus of the present work. Mackey| 



1 



et al. ( 2003k estimated the PISN rate based on t 



leir cal- 
culations of the Pop III star formation rate, predicting 
^2 X 10^ PISNe yr ~ ^ ove r the whole sky above z = 15. 
Weinmann & Lilly (2005) performed a similar analysis 
with more conservative estimates for the star formation 
rate, finding a PISN rate of ^4 yr~^ deg~^ above z = 15 
and ~0.2 yr~^ deg~^ above z = 25, as well as concluding 
that PISNe should be observable out to z = 50 with the 
JWST. 



Subsequent work by Wise & Abel (20051 determined 
the PISN rate based on the collapse of gas into dark mat- 
ter minihalos. Accounting for radiative feedback, they 
concluded that ^0.34 PISNe yr^^ deg~^ were to be ex- 
pected above z = 10, as well as briefly considering the 
dete ctability of PISNe at high redshifts based on models 
from Heger & Woosley (2002). Scannapieco et al. (20051 
presented a more thorough analysis of the visibility of 
PISNe based on a suite of numerical simulations spanning 



hydrodynamics code KEPLER. Mesinger et al. (2006b I 



presented a similar but more general halo mass function 
based analysis, considering the rates and detectability of 
all SNe; they briefly consider primordial PISNe, but fo- 
cus on core-collapse SNe a nd their detectability. More re- 



cently, Trenti et al. (|200^ explored the observable PISN 



rate within the context of the metal-free gas supply dur- 
ing the epoch of reionization. 

Our work improves upon these inv estigations by incor- 
porating updated PISN models from iKasen et al.] 



and determining their observability using the published 
specifications of the Near Infrared Camera (NIRcam) on 
the JWST|^ In the final stages of the work on this pa- 
per w e have become aware of the study by [Pan et al.| 
(20121, who have performed a similar analysis. This pa- 



per addresses the question of PISN observability in a 
nicely complementary way by employing a different nor- 
malization strategy. Different from our assumption that 
viable Pop III pro genitors can only form in unenriched 
minihalos at z > 6, (Pan et al. (2012) derive the star for- 
mation rates required to produce sufficient photons for 
reionization. They then infer the PISN rate correspond- 
ing to different choices for the IMF. Their analysis is 
thus able to probe the charachter of star formation in 
the dwarf galaxies that are the main drivers of reioniza- 
tion, whereas we focus on the minihalos where Pop III 
stars first begin to form. 

This paper is organized as follows. In Section 2 we 
describe our semi-analytic model for the PISN rate. We 
consider the ability of the JWST to detect PISNe at high 
rcdshift in Section 3, and our conclusions are gathered in 
Section 4. Throughout this paper we adopt a ACDM 
model of hierarchical structure formation, using cosmo- 
logica l parameters consisten t with the WMAP 5-year re- 
sults ( [Komatsu et al.||2009[ ): = 0.258; f^A 0.742; 
Vtb = 0.0441; h = 0.719; n/= 0.96; erg = 0.796. 



2. THE PISN RATE 
PISNe are produced only by very massive star s 



(Bromm et al. 



2001 Schaerer 2002 Heger et al 



which are now expected to be rare even for Pop 



2003) 
11 star 



formation. After the first massive star forms, the result 
ing heating from photoionization quickly suppresses the 
density of the remainin g gas in the minihalo, effectively 



density ot tne remainin g gas m the mmihalo, e t tectively 
halting star formation ( "Kitayama et al.' 2004; Whalen 
et al.||20d4| Alvarez et al.|,2006J . The energy released by 
li 



the first PISN dispe rses the gas in the halo and contami 



nates it with metals (iBromm et al. '20(331 Greif et al.|2007 
Whalen et al. 2008; Wise & Abel 200grrGreif et ahliJolOl) 

oi 



Subsequent episodes of star formation are thus delayed 
until the gas is able to re condense into more massive 



cosmological halos (Yoshida et al. 2004 2007 Johnson 



et al.||2007{ [Alvarez et al.„2009| ). While star formation 



will resume at this point, the gas in these systems is ex- 
pected to be enriched beyond the critical metallicity for 
the transition to Population II (Pop II) star formation 
( [Wise fc Ab el"2007', "20081 preif et al.|| 20071 [20081 [2010 ). 

As a result, the stars that form will no longer be massive 
enough to reliably produce PISNe. These explosions are 
thus only expected to occur in minihalos containing pris- 
tine gas that has just crossed the density threshold for 

^ http:/ /www.stsci.edu/jwst/instruments/nircam 
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star formation via H2 cooling, and only one PISN occurs 
per halo. 

It is possible that the first stars formed in binaries or 
small multiples; however, the number of massive stars 



formed per minihalo i s still of order un ity (Stacy et al. 
MolpaTket al .20"TT|[Greif et al.l2011[ ). Given tins, and 
assuming that the time required lor the progenitor star 
to form, live and die is negligible ( Heger et al.||2003 ), we 
can use the formation rate of minihalos to place a robust 
upper limit on the PISN rate. 

The introduction of cosmic feedback has the potential 
to significantly alter this picture. Chemical enrichment 
induces a transition to lower mass Pop II star fomation, 
and thus always has a negative effect on the PISN rate. 
Radiative feedback — especially Il2-dissociating Lyman- 
Werner (LW) feedback — has a more complicated effect. 
The destruction of molecular hydrogen by LW photons 
removes the ability of pristine gas to cool effectively. This 
suppresses star formation, and hence negatively affects 
the PISN rate. However this merely delays star forma- 
tion; as the halo mass continues to increase, the gas even- 
tually becomes self-shielding and proceeds with cooling 
and collapse. The delay effectively increases the amount 
of gas available when star formation begins. This in- 
creases the possibility of multiple massive stars forming 
simultaneously, positively affecting the PISN rate. To 
account for these possibilities, we consider three scenar- 
ios. First we estimate the PISN rate assuming that the 
first stars form unimpeded by cosmic feedback. Then we 
consider a conservative scenario in which both chemical 
and LW feedback affect the PISN rate, but only one star 
forms per halo. Finally, we allow for enhanced massive 
star formation in the LW-affccted halos, and calculate 
the observable rate for each scenario. 



2.1. No-Feedback Limit 

In order to determine an upper limit to the PISN rate, 
we assume exactly one PISN per minihalo, forming as 
soon as the virial temperature of the minihalo exceeds 
the minimum value Td-a required for gas to cool and 
collapse to high densities. We set this to 2200 K based 



on the results of simulations (see § |2.2.1 for details). The 
corresponding critical mass for collapse is given by 



Merit = 10^ M, 







103 K 



3/2 



10 



-3/2 



(1) 



where we assume a mean molecular weight of /i = 1.22, 
appropriate fo r the almost completely neutral IGM at 
high redshifts (Barkana & Loeb 2001). 

We use the analytic Press-Schechter (PS) fo rmalism 
for structure formation ( Press fc Schechter|1974 1 to esti- 
mate the number density Ups of minihalos of mass Af at 
redshift z, given by 



,(M, z) 



2 Pm 
TT M 



din aoiM) 



din M 



./2 



(2) 



where pm is the background matter density, az{M) is 
the standard deviation of overdensities 6 of mass M at 
redshift z, and = 5c/(Jz{M), where Sc is the critical 
overdensity for collapse; the value used here is 6c = 1.686. 



Converting from redshift to cosmic time t(z), where 



tiz) 



1 



dz' 



{i + z')^n^{i + z'Y + n^'' 



(3) 



a first order estimate for the PISN formation rate as a 
function of redshift is given by tt-pisn = dups/dt. How- 
ever, this estimate is only valid while the rate of destruc- 
tion of minihalos via mergers remains small compared 
to the formation rate. The Press-Schechter formalism 
only gives the total number of halos of a given mass at a 
given redshift. As a result, once mergers become impor- 
tant the rate of change of the total number of minihalos 
no longer traces the formation rate. To correct for this, 
we use the expression for the formation rate derived by 
Sasaki] ( (1994 fl 



n+(-) = §«ps(AfcHt,z)^2(§)^: 



(4) 



where D{z) is the growth factor. The PISN rate in this 
upper limit of no feedback, shown in Figure [2] (blue line), 
is then simply given by the halo formation rate: 



fipisNiz) = n+{z). 



2.2. Feedback 



(5) 



The preceding analysis has only nominally incorpo- 
rated the baryonic physics involved through the critical 
mass for H2 cooling. Gas will not successfully cool and 
collapse in all m inihalos that reach the critical mass (e.g., 
Yoshida et al.| [2003). The various feedback mechanisms 
responsible tor this include photoheating from stars in 
nearby halos and the buildup of a background of H2 dis- 
sociating LW photons. Chemical feedback will enrich the 
gas with metals, improving its ability to cool. However, 
gas that is enriched forms lower-mass Pop II stars, effec- 
tively reducing the PISN rate. These feedback mecha- 
nisms can be represented with distinct efficiency factors 
r]{z), such that the true PISN rate will be given by 

ripisfiiz) = r]chcmiz) 77rad(z) n+{z). (6) 

We must include these effects in order to derive a realistic 
estimate for the PISN rate, which we henceforth refer to 
as the conservative feedback case. 

2.2.1. Lyman- Werner Feedback 

LW feedback is of particular importance in any dis- 
cussion of feedback on the first stars as it dissociates 
the H2 molecules primarily responsible for cooling pri- 
mordial gas. This significantly reduces the ability of the 
gas t o cool and works to suppress furth er star forma- 
tion dHaiman et al |pJ97l [Omukai fc Nishi 1999; CiardU 



et al. 



\'}m\ IHaima n et al.||:^00o| |GI3^^F1^ Brand 200l 



2001 




2003 





^002ab; Yoshida et al.| |2003[ lOmukai fc Yoshn 
Mesinger et al. 2006a|~^ As the first stars form 



fluctuations (Barkana & Loeb 2001), we focus here on 
a similarly overdense environment in order to provide a 
conservative estimate of the effects of LW feedback on the 



bee |Mitra et al] l |201l| | for a discussion of the validity of this 
expression. 
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Fig. 1. — a) the factor fux by which LW feedback increases the 
critical mass for star formation; points show the results of the sim- 
ulation, the red line our fit. b) The critical mass for star formation 
in both simulations. Blue crosses mark the critical mass for star 
formation in the absence of LW feedback, red crosses in its pres- 
ence. The blue line shows the best fit critical mass from Equation [l] 
given by a temperature of 2200 K; the red line the critical mass 
from Equation [9] for LW feedback using the fit for /lw given in 
Equation |8] The green line marks the halo mass corresponding to 
a virial temperature of 10* K, and the black dashed line represents 
the critical mass employed, accounting for atomic cooling in halos 
with virial temperatures above 10* K. 



PISN rate. While representative of the LW background 
in Pop III star formation sites, this is likely not repre- 
sentative of the 'avera ge' LW background in . the universe 
(e.g., |Machace k et aLpOOlj [Mesinger et al.|2006a^ . 

We will describe our radiative feedback simulations in 
detail elsewhere, and present only a brief summary here. 
We employ a set of two cosmological simulations using 
a modified version of the A^-Body/TreePM S PH code 
GADGET ( |Springel|[2005l [Springel etaLl[200l] ) to gauge 
the effects of LW feedback on the PlSjSi rate. These 
simulations, carried out in a box of size 3.125 co- 
moving Mpc and starting from a redshift of z = 127, 
were performed to investigate the formation of the first 
dwarf galaxies in 10^ M© halos at z = 10. In order to 
obtain high resolution in the halo containing the first 
galaxy while retaining information about structure on 
large scales, a zoomed simulation technique was used, 
with the highest resolution dark matter (gas) particles 
having a mass of 2350 (484) Mq . This allows halos with 
masses > 2 x 10''^ M©, to be resolved with > 100 dark 
matter particles. 

The first of the simula tions we employ is si milar to sim- 



ulation Z4 presented in Pawlik et al. (2011 1, and follows 



the non-equilibrium chemistry and cooling of the primor 
dial atomic and molecular gas, including star formation 
but not the associated feedback. It thus provides a useful 
reference against which the effects of LW feedback can 
be discussed. Here, gas particles were turned stochasti- 
cally into star particles at densities nn > 500 cm^'^ on 
a dynamical time scale. Star particles were considered 
simple stellar populatio ns using the zer o-metallicity top- 
heavy IMF models from Schaerer ( 2003 ) . Henceforth this 
simulation is referred to as Simulation A. 
The second simulation employed here, which we re- 



fer to as Simulation B, is identical to Simulation A ex- 
cept for the inclusion of LW feedback. Galculation of the 
feedback was carried out by considering the contribution 
from both star particles and a uniform LW background. 
The combined LW background is normalized to approx- 
imate the LW background evolution shown in |Greif fc| 



Bromm ( 2006 1 in the optically thin limi t, but with the 



application of a self-shielding correction (Wolcott- Green 
etaL][20TT|) . 

'i'he etticiency of LW feedback, yyLw, can be expressed as 
the ratio of the formation rate of minihalos at the critical 
mass with LW feedback to that without: 



n+(Mcrit,LW,-Z) 

n+{Mcrit,z) 



(7) 



This requires determining the factor /lw(-z) by which LW 
feedback increases the critical mass for star formation. 
We do this by determining the critical mass required 
for stars to form in both simulations, as traced by the 
lowest mass halo that is actively forming stars for the 
first time. The resulting critical mass in each case is 
shown in Figure [Tjs; blue points denote the critical mass 
in Simulation A, red points in Simulation B. Without 
any negative feedback, we expect stars to form when the 
halo reaches the critical mass given in Equation [T] We 
find the that the resulting critical mass in Simulation A 
is best fit by a critical virial temperature of 2200 K, as 
shown in Figure [T]d. Note that this includes the effects of 
dynamical heating; gas in isolated halos would collapse 
and form stars at lower masses. 

Determining the critical mass for star formation in the 
LW feedback simulation in the same manner as above, 
we can then determine /lw, shown in Figure [T^. We find 
that /lw is well fit by a functional form of 



f^^.{z) = -6.23 X 10"^ erf[0.094(z - 0.204)] 
-1-6.23 X 10^ 



(8) 



for z > 10, where erf(z) is the error function. Afcrit.Lw is 
then given by 



Mr, 



/lw Mc, 



(9) 



shown in Figure [Tja. When the virial temperature of the 
halo reaches 10** K, cooling via atomic hydrogen becomes 
efficient and molecular hydrogen becomes unimportant. 
Once this threshold is reached, LW feedback cannot fur- 
ther suppress star formation, and Merit, lw holds steady at 
a constant virial temperature of 10^ K. This is refiected 
as an increase in the formation rate of halos affected only 
by LW feedback, as the formation rate of atomic cooling 
halos is increasing prior to z = 10. This transition can 
be clearly seen in Figure [2] marked by a sudden jump 
in the halo formation rate. As ionizing photons cannot 
easily escape the immediate vicinity of the star that pro- 
duced them, their impact on neighboring halos is small 
compared to that of LW photons at the highest redshifts. 
We thus set Tyrad — ?7lw for simplicity. This approxima- 
tion becomes increasingly unphysical close to the epoch 
of reionization at z ~ 10. 

2.2.2. Chemical Feedback 

The process of chemical enrichment is another crucial 
factor for determining the PISN rate. Gas that has been 
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Fig. 2. — a) npisN in the upper limit of no feedback (blue), with 
chemical feedback (green), LW feedback (red) and the resulting 
PISN rate for the conservative (chemical plus LW) feedback case 
(black), b) Same as (a), but for enhanced massive star formation. 



enriched beyond a critical metallici ty of Zait ^10 
will no longer fo rm Pop III stars (pro nirn et al. 



2001 



Schneider et al.||2002| fBromm fc Loeb||2U03| ), and hence 
no PibiNc. Chemical teedback can thus be represented 
as the fraction of halos forming from pristine gas at a 
given redshift. Realistic three-dimensional simulations of 
this process starting from cosmological initial conditions 
have become possible in the past decade, showing that 
enrichment by Pop III SNe, if they are highly energetic, 
proceeds very inhomogeneously, enriching the IG M be- 
fore penetrating into denser regions (Scannapicco et al 



2005] [( 

^ra]|2008; M aio eraTpOTOf 
In modeli ng rychem, we use the results of IFurlanetto fc 



Greif et al.|[2007l |Tornatore et al.„2007i |Wise 



|Loe b| ([2005 1 . Their semi-analytic treatment of SN winds 
utilizes the Sedov (1959) solution for an explosion ex- 



panding into a uniform medium and yields a probability 
function Ppristine(^) that the gas in a newly formed halo 
is pristine. This is plotted in Figure 2 of their paper for 
various strengths of chemical feedback. Wc identify this 
quantity as the fraction of newly collapsed halos that 
have been polluted with metals, ??chom- Given the re- 
cent detection of pristine gas at z = 3 by | Fumagalli| 
et al. (2011), we choose the weakest fe edback scenario 
presented by Furlanetto & Loeb (2005) among the sce- 
narios that incorporate a clustering ot sources. The re- 
sulting PISN rate is given by the green line in Figure [2j 

2.3. Enhanced Massive Star Formation 

Gas cooling and subsequent star formation in halos af- 
fected by LW feedback can be delayed until nearly an or- 
der of magnitude more gas is available for star formation 
(Figure [TJ. This increases the likelihood that multiple 
massive stars form per halo, offseting the negative effects 
of LW radiation considered above. We quantify this by 
positing that the number of PISNe produced per halo 
at redshift z is given by the ratio of the critical mass in 



10' 



10' 



■g 10-1 



Miralda-Escude & Rees 97 



' Mackey et al. 03 



Weinmann & Lilly 05 
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Fig. 3. — The observable PISN rates in number per year per 
JWST field of view above a given redshift in the upper limit of no 
feedback (blue line), in the conservative feedback case (solid red 
line), and the enha nced star formation case (dashed red line). The 
rates calculated by 'Miralda-Escude & Rccs ( 1997^, iMackey et al.| 
(2003), Weinmann & Lilly (2UU5) and Wise & Abel (2UU5) are also 
shown for reference. Ked points account tor teedback; blue points 
do not. 



the presence of LW feedback Merit, lw to the critical mass 
in the no- feedback case Merit- For example, at z = 17, 
Merit, Lw/-^crit ~ 1-4, SO for evcry 10 pristine halos that 
form, 14 PISNe are produced. In this case the PISN rate 
is modified such that 



/ X Merit, Lw(2;) / N r \ ■ r \ 

Merit (^) 



(10) 



The resulting enhanced PISN rate can be seen in 
Figure (2)3. In contrast to the conservative feedback case, 
the net effect of LW feedback is much less significant here, 
with chemical feedback controlling the final PISN rate. 

2.4. The Observable Rate 

The observed PISN rate per unit time per unit redshift 
per unit solid angle is given by 



dN 
dtohs dz dft 



dN dV 
dtohs dV dzdfl 



1 



dN 



, dr 



(11) 



(1 + z) dtr 



.dV ^ dz' 



Cosmological time dilation between iobs a-nd t^m is ac- 
counted for by the (1 -I- z) in the denominator; dV is 
the comoving volume element and r(z) is the comoving 
distance to redshift z given by 



r{z) 



dz' 



Ho Jo ^n„,{i + z')^ + n^' 



(12) 



where c/Ho is the Hubble distance. With the assump- 
tions outlined above, we estimate the PISN rate in events 
per year per comoving Mpc'^ in the source rest frame: 



dN 



dtrrr, dV 



'^pisn(^)- 



(13) 



These results — shown in Figure |3] — are in reasonable 
agreement with previous work; our no-feedback limit of 
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one PISN per minihal o is somewhat mor e cons ervative 
than that employed by Weinmann & Lilly (2005), but in 
general agreement. Likewise, our conservative feedback 
rate is in go od agreement with the rate found by [Wise "fc] 
Abel (2005), which also accounted for feedback. The dis- 



crepancy wi th the remaining rat es can be attributed to 
the fact that Mackey et al. ( 2003 1 employed an optimistic 
star formation ehicie ncy ot 77* = 0.10, while |Miralda" 



Escude & Rees ( 1997) performed a rough estimate of the 
all-sky supernova rate based on the observed metallicity 
of the IGM. This estimate of course includes Population I 
and II supernovae in addition to PISNe, and is included 
only for reference. 

3. JWST OBSERVABILITY 

While PISN explosions are predicted to be extremely 
energetic, the highest redshift events will still be unob- 
servable, and those at lower redshifts will be above the 
detection limits of the JWST for only a fraction of their 
lifetimes. To determine the observability of these explo- 
sions we must consider both how bright they will be at 
a given redshift and how long they will remain visible. 
To span the uncertainties arising from variation in the 
progenitors of PISNe we consider a set of four models 
from Kasen et al.i ( |2011| ), namely their R250, B200, R175 
and HelOU models. R250 and R175 are red supergiants 
of 250 M0 and 175 Mq, respectively, spanning the mass 
range of succesful explosions. We also consider a more 
compact 200 Mq blue supergiant (B200) and a 100 M© 
bare helium core (HelOO). All models considered here die 
as PISNe, with explosion energies ranging from 7 x 10^^ 
ergs (R250) to 2 x 10^^ ergs (R175). The late-time lumi- 
nosity is powered by the decay of ^^Ni produced during 
the explosion. 40 Mq of ^^Ni are produced in the R250 
model, while only 5, 2 and 0.7 M© are produced in HelOO, 
B200 and R175, respectively. 

Bare helium cores, such as the HelOO model, arc of 
particular interest in light of recent work finding that 
rapidly rotating stars can encounter the pair-production 
instability in progenitors with masses as low as 65 M© 
( [Chatzopoulos fc Wheeler|2012||Yoon et al.|2012] . Com- 
bined with recent work iinding that Pop 111 stars are 
both less mass ive and more rapidly rotating than previ 



ously thought ([Stacy et al.|2010l[20]H|2012|park et al. 
20T1I [Greit et al.||2(jll[ |20J2[ ), tte explosion ot a rapidly 



rotating helium core formed by homogeneous evolution 
represents an intriguing possibility. 

We first describe our technique for fitting a blackbody 
to these four PISN models before considering their visi- 
bility with the JWST. We then estimate the total observ- 
able number in each case and for each feedback prescrip- 
tion. Finally, we briefly discuss the challenges involved 
in actually identifying PISNe as such. 



3.1. A Simple Lightcurve Model 

In order to determine how long a PISN will be visi- 
ble we must model the source spectrum. Given the large 
mass involved, the ejecta will remain optically thick until 
late times, so we make the reasonable assumption that 
the PISN emits as a blackbody for the majority of its visi- 
ble hfetime. Using the U, B, V, R, I, J, H, and K absolut e 
magnitude light curves presented inlKasen et al. (2011), 
we perform a least-squares fit to find the combination 




100 150 200 250 300 350 400 450 
Time 



Fig. 4. — The temperature and radius of our blackbody fits as a 
function of source frame time since the explosion for models R250, 
B200, R175 and HelOO. The secondary rise in temperature seen in 
R250 and B200 is caused by the decay of ^^Ni reheating the ejecta 
at late times. Note also that the apparent radii begin to decrease 
again at late times; this can be interpreted as the photosphere 
receding into the ejecta as the material cools. 



of temperature T and radius i?sN that best matches the 
broadband magnitudes at each point in time. This is 
done with the assumption that the specific luminosity of 
the PISN LgN,A at wavelength A is given by 



(14) 



where Bx is the Planck function. The resulting fits for 
the evolution of the temperature and radius of the PISN 
models are shown in Figure |4] Note that our blackbody 
assumption breaks down at late times when the photo- 
sphere begins to recede into the ejecta. This is mani- 
fested as an apparent decrease in the radius of the PISN 
remnant. 

With this information we can then calculate the spe- 
cific flux -Fx. cm iu the rest frame and — accounting for red- 



shift and cosmological dimming- 
for a source at redshift z: 



-in the observer's frame 



Fx ,obs 



Rrk 



1 2 



Bx'jT) 
1 + z ' 



(15) 



Here A' — A/(l -I- z) accounts for redshifting, and the lu- 
minosity distance = (1 -I- z)r{z) accounts for cosmo- 
logical dimming. Convolving this spectrum with a fllter 
function (j>x{^) yields the observable flux in filter X: 



-Fobs.: 



0x(A)FA,obs(r, -RsN, z)dX. 



(16) 
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Fig. 5. — Left: Lightcurves for the |Kasen et al.| | |2011| l R250 (top) and B200 (bottom) models as they would be observed by JWST's 
F444W NIRCam filter at 2 = 5, 10, 15, 20, 25 and 30. The flux limits for a lO'' s (dashed line) and lO" s (dotted line) exposure are shown 
for reference. Right: The visibility time Atvis in years for R250 (top) and B200 (bottom) as a function of redshift for each of the NIRcam 
wide filters. Note that the axes are scaled independently. Similar plots for models IfelOO and R175 are included in the appendix. 



3.2. Visibility 

The NIRCam instrument on the JWST will observe 
the early universe through a number of narrow, medium- 
width, and wide filter^ The widest, longest- wavelength 
filter, F444W, will observe from 3.3 to 5.6 ^m with a 
sensitivity limi t of 24.5 nJv required for a lOcr detection 
in 10^ seconds ( Gardner^ et al. 2006). Shown in the left- 
hand column of l^'igure JSlis the observable flux as it would 
appear in the F444W NIRCam filter at various redshifts 
for the most and least easily observable models, R250 
and B200, respectively. See Figure [7] for why these two 
were chosen; models HelOO and R175 can be found in the 
appendix. The flux limits for the filter of 4.4 x 10~^^ erg 
s~^ cm~^ for a 10^ s exposure and 4.4 x 10""'^^ erg s~^ 
cm~^ for a 10** s exposure are also shown for reference. 
We see that the brightest explosions (R250) would be 
visible to beyond z ^ 25, but are never so bright as to 
be detectable with current generation telescopes. This is 



consistent with the non-detection by Frost et al.| ( |2009 1 
in a search of the Spitzer/IRAC Dark J^ield tor possible 
Pop III PISN candidates. 

To account for absorption of flux by neutral hydrogen 
along the line of sight we implement a simple model of 

^ http: / / www.stsci.odu/jwst / instruments / nircam / instrument- 
design/filters 



instant reionization at z = 10. For sources above this 
redshift, we assume no flux is observed shortward of the 
rest frame Lya line. This is not relevant for the F444W 
NIRcam filter as Lya does not redshift into the filter until 
z 40, when the lightcurve is already far below even the 
10^ s sensitivity limit. It does however have an effect, 
albeit a smafl one, on the F115W and F090W filters. 

At low reds hifts the dur ation of the lightcurve pre- 
sented in Kascn et al. (2011 ) is not quite long enough for 
the observed flux to reach the sensitivity limit; we ex- 
tend it to the limit by extrapolating assuming a power- 
law scaling. The visible time Atvis is then simply given 
by the time the lightcurve is above the filter sensitivity 
limit. Shown in Figure [5] are the visibility times as a 
function of redshift for each of the NIRcam filters. 

3.3. The Observable Number 

With this estimate for Afvisi we may finally calculate 
the observable number of PISNe on the sky, given by the 
product of the PISN rate at z, as seen in the observer 
frame, and the time a PISN at z is visible, A^vis- This 
yields an estimate for the number of PISNe visible on 
the sky at any given time per unit redshift per unit solid 
angle: 

dN ^ dN 
dz dil (itobs dz dil 



8 HUMMEL ET AL. 



R250 B200 




Fig. 6. — Upper and lower limits for the number of PISNe per JWST FoV above redshift z with different feedback prescriptions. The 
observable numbers for a 10® s exposure assuming the R250 model are shown on the left; the B200 model is employed on the right. Solid 
blue lines show an upper limit to the observable number in the case of no feedback, solid red lines an estimate for the observable number 
in our conservative feedback scenario, and dashed red lines the number in the enhanced star formation case. Note that the x-axis is scaled 
independently in each panel. Also shown in the R250 panel are the observable numbers for a lO** s exposure for the no-feedback (dot-dashed 
blue), conservative feedback (dot-dashed red), and enhanced star formation (dotted red) scenarios. The B200 model is not visible above 
z = 5 in a 10* s exposure. 



Figure [6] shows the number of PISNe per JWST field of 
view (FoV) above redshift z in a 10^ s exposure for all 
three feedback cases for models R250 and B200. Mod- 
els HelOO and R175 can be found in the appendix. For 
the R250 model the results for a 10'* s exposure are also 
included; B200 is not visible at high redshifts in a 10^ s 
exposure. 

In the optimistic case of an R250-type PISN with no 
feedback we expect -0.2 PISNe per JWST FoV for a 
10® s exposure. In the most pessimistic case of a B200- 
type PISN with strong negative feedback this number 
drops to ~2.5xl0~^ per FoV. The actual number de- 
tected by the JWST will most likely lie somewhere within 
this range. Given this, we conclude that a single deep 
pencil-beam survey is unadvisable for detecting PISNe, 
as there aren't enough in a given field to ensure a detec- 
tion, even in the most optimistic upper limit. This sug- 
gests that a mosaic search, covering a larger area with 
shorter exposure times, may be the best approach to en- 
sure finding a Pop III PISN. 

3.4. PISN Identification 

The exceedingly long duration of their lightcurves 
poses a serious challenge for identifying PISNe. When 
combined with the cosmological time dilation factors in- 
volved, PISN lightcurves can last for decades; the highest 
redshift events will last longer than the projected mission 
lifetime for the JWST, making the detection of PISNe 
by searching for transients difficult at best. However, 
a multi-year campaign might be able to detect photo- 
metric variations; for example, the HelOO model would 
appear to decline in bright ness by ^0.3 m agnitudes per 
year at redshift z ~ 10 (Kas en et al.|2011 ). Additionally, 
PISN colors become redder over time as the photosphere 
receeds into the ejecta and metal line blanketing sup- 
presses flux in the bluer bands. While likely insufficient 
to unambiguously identify PISNe, this could be useful in 
selecting candidates for spectroscopic follow-up. 

While the peak bolometric luminosities and spectra of 
PISNe resemble those of typical Type la and Type II su- 



pernovae ( [Joggerst k Whalen||201l| [Kasen et al^|2011|). 
relatively l ittle mixing occurs during tlie explosionl J og- 
gerst & Whalcn 20Tl| iChen et al. 1201 1|). As a result 



PlSiNe mostly retain their omon-iayer structure during 
the explosion, and metal lines do not appear in the spec- 
trum until late times when the photosphere has receeded 
deep into the ejecta. Some lighter elements may appear, 
but the early sp ectrum of a PISN will be d evoid of Si, 
Ni and Fe lines ("Joggerst & Whalen 2011 1. This may 
provide the spectroscopic signature needed to identify 
PISNe. 

4. DISCUSSION AND CONCLUSIONS 

In this work, we have examined the source density of 
PISNe from Pop HI stars and considered their detectabil- 
ity with the JWST. We conclude that the limiting factor 
in detecting PISNe will be the scarcity of sources rather 
than their faintness , in agr eement with the conclusions of 
Weinmann fc Lilly] p005| . The brightest PISNe should 
be readily detectable with the longest wavelength NIR- 
cam filters out to z ~ 25; the problem is the overall 
scarcity of sources. 

We have derived an estimate for the observable PISN 
rate, finding an upper limit of just over 0.01 PISNe per 
year per JWST FoV in the case of negligible chemical 
and radiative feedback. We also find that the inclusion 
of feedback can reduce the PISN rate by an order of mag- 
nitude, to ^0.001 per FoV. Accounting for the possibil- 
ity of enhanced massive star formation in halos affected 
by LW radiation improves this rate slightly, to ^0.003 
per FoV. The derived PISN rate then allows us to place 
an upper limit on the observable number of PISNe in a 
10® s exposure of -0.2 PISNe per JWST FoV in the no- 
feedback case. The most pessimistic case of a B200-type 
PISN with strong negative feedback reduces this number 
to -2.5 X 10--* per FoV, or one PISN per 4000 JWST 
fields of view. 

The long duration of PISN lightcurves imply that spec- 
troscopic follow-up of PISNe will likely be of great im- 
portance. PISN lightcurves can last for decades when 
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Fig. 7. — The observability of the R250, HelOO, B200, and R175 models (clockwise from upper left) from |Kasen et al.| l [20lT|l using the 
JWST's NIRcam F444W filter. Shown is the possible range for the number of JWST FoVs required to detect iU sources as a function of 
exposure time. The blue range is for all PISNe, and the red for PISNc from z > 15. The lower boundaries correspond to the no-feedback 
upper limit to the PISN rate and the upper boundaries to the conservative feedback rate. From left to right, the black lines represent the 
number of pointings possible in a total of 10'', 10'', and 10* s for a given exposure time. 



combined with the cosmological time dilation factors at 
high redshifts, making the detection of PISNe by looking 
for transients untenable. However, a multi-year cam- 
paign could identify candidates photometrically, and the 
lack of metal lines in the spectrum at early times could 
provide a spectroscopic signature for identification. 

We find that the main obstacle to observing PISNe is 
the paucity of sources. Beyond a moderate exposure time 
of a few times 10"* s, the observability of bright PISNe is 
not a strong function of exposure time and is instead 
controlled by the source density; this is evident from 
Figure [T] where we have shown the number of JWST 
FoVs required to detect 10 PISNe (blue) as a function 
of exposure time for each of our PISN models. The up- 
per boundaries correspond to our conservative estimate 
for the PISN rate in the presence of feedback, and the 
lower boundaries to the estimated rate without feedback. 
We can see that the decrease in the required number 
of JWST pointings slows considerably beyond ^lO'^ s, 
hence a deep pencil-beam survey would not be advis- 
able in searching for PISNe. Even for only high-redshift 
sources (z > 15; red) the dependence on exposure time 
is still minimal, being controlled by the lack of sources 
once the required imaging depth is reached. 

Of particular interest in Figure [7] are the black lines 



representing the total number of pointings possible in 
10®, 10*", and 10* s for a given exposure time and their 
location relative to the observability ranges in blue and 
red. 10® s is approximately the limit of what would be 
possible with a dedicated deep-field campaign; lO'^^s is 
the limit of the observations the JWST could make in a 
year assuming NIRCam is in use one third of the time; 
10* s (~10 years) is the projected mission lifetime. 

While the detection of a PISN from a 'first' star at very 
high redshifts would be exciting and is in fact possible 
given the detection limits of the JWST, the scarcity of 
sources at these redshifts means that such a detection 
would be highly contingent on serendipity. Even in the 
most optimistic case, with all available minihalos pro- 
ducing an R250-type PISN, the observability range for 
such events lies well above what is possible even in a full 
year of observations, though a few may be detected over 
the lifetime of the telescope. The detection of a PISN 
at lower redshifts appears to be more realistic. As the 
faintest PISNe (R175 and B200) are effectively unobserv- 
able, PISN searches should focus on looking for PISNe 
similar to the R250 and HelOO models. In this case, 
the strategy with the highest likelihood of detection will 
be a mosaic survey of many moderately deep exposures. 
This is clear from Figure [8j where we show the number 
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JWST Fields of View JWST Fields of View JWST Fields of View 

Fig. 8. — The total number of PISNe observable with a campaign of 10®, lO'^ and 10* s (from left to right) as a function of survey area for 
the R250 PISN model. In each case, the total campaign time is apportioned equally over the total survey area to determine the exposure 
time for individual pointings. The blue region represents all PISNe, the red only PISNe from z > 15. Upper boundaries correspond to the 
no-feedback upper limit to the PISN rate and lower boundaries to the conservative feedback case. For reference we mark the case of only 
one PISN visilalc (dashed line). 



of PISNe that will be observable with the JWST in ob- 
serving campaigns totalling 10^, 10'' and 10^ s for the 
R250 PISN model. The exposure time for each point- 
ing varies with the total area covered by the survey in 
order to keep the total observing time constant. Up- 
per boundaries correspond to the number visible in the 
no-feedback case, lower boundaries to the conservative 
feedback case. As in Figure [7j the blue region shows the 
observable number from all redshifts, the red region only 
those from z > 15. We see that the observable number 
increases until the resulting exposure time is no longer 
sufficient to detect PISNe. The optimal search strategy 
then will be to cover as large an area as possible, going 



only as deep as necessary, possibly in a similar manner to 
the ongoing Brightest of Rei onizing Galaxies survey with 



the Hubble Space Telescope (Trenti et al. 2011 
eraL][20T2| . 
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APPENDIX 



Included here are the observable properties of the HelOO and R175 models from Kasen et al. (20111, including the 
observable flux and time visible (Figure [lO]) and the observable number for each feedback scenario (I'lgurejO]). These 
plots are to be compared to Figure ^ anoFigure [6] in §[3j Both cases may be detected with NIRCam beyond z = 15 
for deep exposures of 10^ s. 

Also included is an extended version of Figure |8] showing the detectable number in observing campaigns totalling 
10*^, 10^ and lO^s for all four PISN models (Figure [TT]) . Only R250-type PISNe will be detectable above reshift 15, 
and R175- and B200-type PISNe will only be detectable in observing campaigns totalling greater than 10'' s. 
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Fig. 9. — Upper and lower limits for the number of PISNe per JWST FoV above redshift z with different feedback prescriptions. The 
observable numbers for a 10^ s exposure assuming the HelOO model are shown on the left; the R175 model is employed on the right. Solid 
blue lines show an upper limit to the observable number in the case of no feedback, solid red lines an estimate for the observable number 
in our conservative feedback scenario, and dashed red lines the number in the enhanced star formation case. Note that the x-axis is scaled 
independently in each panel. 
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Fig. 10. — Left: Lightcurves for the |Kasen et al.| | [2011| | HclOO (top) and R175 (bottom) models as they would be observed by JWST's 
F444W NIRCam filter at z = 5, 10, 15, 20, 25 and 30. The flux limits for a lO'' s (dashed line) and lO" s (dotted line) exposure are shown for 
reference. Right: The visibility time Atvia in years for HelOO (top) and R175 (bottom) as a function of redshift for each of the NIRcam 
wide filters. Note that the axes are scaled independently. 
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Fig. 11. — The total number of PISNe observable with a campaign of 10^, lO'^ and 10* s (from top to bottom) as a function of total survey 
area for a given PISN model. In each case, the total campaign time is apportioned equally over the survey area to determine the exposure 
time for individual pointings. The blue region represents all PISNe, the red only PISNe from z > 15. Upper boundaries correspond to the 
no-feedback upper limit to the PISN rate and lower boundaries to the conservative feedback case. The dashed line marks one PISN visible. 
Note that not all axes are scaled the same. 



